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ABSTRACT 

Context. Precise analyses of stellar radial velocities reveal intrinsic effects causing wavelength shifts of spectral lines (other than 
Doppler shifts due to radial motion), such as gravitational redshifts and convective blueshifts. 

Aims. Gravitational redshifts in solar-type main-sequence stars are expected to be some 500 m s -1 greater than those in giants. Such 
a signature is searched for between groups of open-cluster stars which share the same average space motion and thus have the same 
average Doppler shift. 

Methods. 144 main-sequence stars and cool giants were observed in the M67 open cluster using the ESO FEROS spectrograph, 
obtaining radial velocities by cross correlation with a spectral template. Binaries and doubtful members were removed, and averaged 
low-noise spectra calculated for different classes of stars. 

Results. M67 dwarf and giant radial-velocity distributions are well represented by Gaussian functions, sharing the same apparent 
average radial velocity within =s 100 m s _1 . In addition, dwarfs in M67 appear to be dynamically hotter (tr = 0.90 km s -1 ) than giants 
(o- = 0.68 km s" 1 ). 

Conclusions. Explanations for the lack of an expected signal are sought: a likely cause is the differential wavelength shifts produced 
by different hydrodynamics in dwarf and giant atmospheres. Radial-velocity differences measured between unblended lines in low- 
noise averaged spectra vary with line-strenth: stronger lines are more blushifted in dwarfs than in giants, apparently 'compensating' 
for the gravitational redshift. 

Synthetic high-resolution spectra are computed from 3-dimensional hydrodynamic model atmospheres for both giants and dwarfs, 
and synthetic wavelength shifts obtained. In agreement with observations, 3D models predict substantially smaller wavelength-shift 
differences than expected from gravitational redshift only. The procedures developed could be used to test 3D models for different 
classes of stars, but will ultimately require high-fidelity spectra for measurements of wavelength shifts in individual spectral lines. 

Key words. Stars: fundamental parameters - Open clusters and associations: individual: M67 - Stars: late-type - Stars: atmospheres 
- Gravitation - Techniques: radial velocities 



1. Introduction 

1.1. Gravitational redshift 

Spectroscopic measurements of stellar radial velocities have 
reached levels of precision which not only enable new types of 
studies, but also require a better understanding of effects other 
than a star's barycentric motion that may cause displacements 
of spectral-line wavelengths. One such effect is the gravitational 
redshift originating from the propagation of light between differ- 
ent gravitational potentials at the source and at the observer. 

Any theory involving conservation of energy requires that 
radiation leaving a gravitational field loses energy, and with the 
early 20 th century understanding of the photon and the relation 
between its energy and wavelength, predictions could be made 
of the wavelength displacement expected for radiation leaving 
various stars, and these effects have been sought observationally. 
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While a determination of an absolute wavelength shift re- 
quires the true stellar center-of-mass motion (and its induced 
Doppler shift) to be known, relative shifts can be determined 
between stars sharing a common space motion. This applies to 
the orbit-averaged velocity of binary stars (where each compo- 
nent shares the same systemic velocity), and for stars in [open] 
clusters which share the same space velocity vector. 

The first claimed confirmation of the predicted gravitational 
redshift came from the measurement of the apparent radial ve- 
locity of Sirius B, the white-dwarf companion to Sirius (Adams 
1925), identified as the difference from the motion expected in 
its 50-year orbit. However, this faint white dwarf is overwhelmed 
by scattered light from the much brighter primary and the valid- 
ity (and even the honesty!) of these early measurements has been 
questioned (Hetherington 1980; Greenstein et al. 1985). 

Gravitational redshifts have since been more reliably de- 
termined for numerous white dwarfs in binary systems and in 
nearby open clusters such as the Hyades and the Pleiades, find- 
ing values typically on order 30-40 km s _1 , as measured against 
wavelength positions in the spectra of ordinary stars. Also, very 
much larger gravitational redshifts have been identified from 
spectra of X-ray bursts on neutron-star surfaces. 
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The gravitational redshift is proportional to the gravitational 
potential at the stellar surface and scales proportional to the stel- 
lar mass, and inversely as its radius R: 

Vgrav = GM/Rc, 

where G is the gravitational constant. The expected solar 
gravitational redshift for light escaping from the solar photo- 
spheric surface to infinity is 636.486 ± 0.024 m s _1 (Lindegren 
& Dravins 2003) and - since the shift diminishes with distance 
from the stellar center as r~ x - becomes 633.5 m s _1 for light in- 
tercepted at the Earth's mean distance from the Sun (r = 215/? ). 
A solar spectral line instead formed at chromospheric heights 
(30 Mm, say; r = 1 .04R Q ) will have this redshift decreased by 
some 20 m s _1 , and a coronal line by perhaps 100 m s _1 . 

For other stars, the shift scales as (M/MqXR/Rq) -1 , or as 
(g/g )(R/R Q ), where g is the surface gravity. For distant stars, 
additional lineshifts may originate from large-scale gravitational 
fields in the Milky Way Galaxy. For instance, spectra of stars 
in the central Galactic bulge (R* 1 kpc) may be gravitation- 
ally redshifted by 300-400 m s _1 , while stars in the Magellanic 
Clouds (R* =! 55 kpc) might be blueshifted by a similar amount, 
as seen by an observer near the solar galactocentric position at 
R b s - 8.5 kpc (Lindegren & Dravins 2003). 

Across the Hertzsprung-Russell diagram, the gravitational 
redshifts are not expected to vary much along the main se- 
quence between A5 V and K0 V (being around 650 m s _1 ), 
but to reach twice that value for more massive early-B stars 
with M 10M o . For giants the shifts decrease substantially 
and for red supergiants it should be only some 20 or 30 m s _1 , 
three orders of magnitude smaller than for white dwarfs (Griffin 
1982; Dravins 2005). Thus, the difference between the shifts in 
dwarf and giant stars is expected to be some 0.5 km s _1 , or- 
ders of magnitude greater than current measurement precision. 
Such values are significant, both to study the precise dynamics in 
binary- and multiple-star systems (Griffin 1982, Pourbaix et al. 
2002), or to disentangle wavelength shifts caused by stellar at- 
mospheric dynamics by convective motions (e.g., Dravins 1982; 
Asplund 2005; Dravins et al. 2005). Together with photospheric 
line asymmetries, such intrinsic wavelength shifts make up sig- 
nificant diagnostic tools for probing the 3-dimensional structure 
of stellar atmospheres (Dravins 2008; Nagendra et al. 2009), a 
field where additional diagnostics are needed to reliably segre- 
gate among hydrodynamic models predicting significantly dif- 
ferent chemical compositions, oscillation properties, and so on. 

1.2. Searching for shifts in ordinary stars 

At first sight, the Sun would appear to be a promising source 
for determining absolute gravitational redshifts, since the radial 
Sun-Earth motion is both very small and well known. However, 
already long ago, spectral lines in the Sun were found to be 
shifted by different but comparable amounts, including varia- 
tion of absorption-line wavelengths across the solar disk from 
its center towards the limb that was initially a mystery. This is 
now well understood as convective lineshifts, originating from 
a summation of contributions from hot and rising (thus locally 
blueshifted) elements, and cool and sinking ones, and various 
effects of observational perspective across the disk. Many ef- 
forts were made to measure the solar gravitational redshift (and 
several publications claimed its detection) although, in hind- 
sight, almost all such studies were constrained by the limited 
understanding of the interplay with other mechanisms causing 
subtle wavelength shifts (besides effects from atmospheric mo- 
tions, also the limited accuracy for laboratory wavelengths, and 
even effects of pressure shifts come to play). Possibly, the most 



convincing demonstration of a gravitational redshift in the so- 
lar spectrum was by Beckers (1977). The magnetic fields in 
the cores of sunspots inhibit convective motions and greatly re- 
duce the amount of both convective wavelength shifts, and their 
center-to-limb variation, enabling rather precise wavelength shift 
measurements in magnetically insensitive lines. From such lines, 
Beckers deduced a shift of 613 m s _1 , consistent with the ex- 
pected magnitude. 

Just as for white dwarfs, systematic wavelength displace- 
ments between ordinary stars may be identified in binaries shar- 
ing a common systemic velocity, and for stars in open clusters, 
moving in space with sufficiently small velocity dispersion. The 
earliest such attempts seem to have been by E. Finlay-Freundlich 
who, starting around 1914, tried to find statistical evidence of 
gravitational redshifts (for a lucid review, see Hentschel 1994). 
Somewhat later, Trumpler (1935), recognized that massive O 
stars should be excellent targets, because they have the largest 
gravitational redshift while they often appeared in young clus- 
ters where they, being massive, by equipartition of energy should 
have the smallest velocity dispersion. In some clusters, shifts of 
some 10 km s _1 were found relative to comparison stars, sug- 
gesting enormous stellar masses. A later understanding of these 
stars showed that shifts should be much less (Conti et al. 1977), 
and while effects of stellar variability, line asymmetries, and 
differences of apparent velocities in different spectral features 
did begin to get appreciated, also a certain literature grew on 
suggested 'anomalous' redshifts in binary high-luminosity stars, 
and other objects. 

Only after measuring precisions for individual stars had 
reached levels of around 100 m s _1 , did it become practical to 
search for systematic differences in apparent velocities among 
ordinary stars, e.g., between main-sequence dwarfs and giants. 
One hint came from a study of the open cluster NGC 3680 by 
Nordstrom et al. (1997), who found the giants (though the sam- 
ple contained only six such stars) to be blueshifted on average 
0.4 km s _1 relative to the dwarfs, as expected from gravity. 

For certain groups of stars, including those in undispersed 
moving clusters, absolute stellar radial motion can be deter- 
mined from precise astrometry, not invoking any spectroscopy. 
The difference to spectroscopic lineshifts then reveals shifts in- 
trinsic to stellar atmospheres. In the Hyades (Madsen et al. 
2002), such signatures were seen for some classes of hotter 
main-sequence stars, but the spectral displacements of the (only 
three) Hyades K-giants did not show any significant deviation 
from that of main-sequence K-type stars. 

Different gravitational redshifts in dwarf and giant stars have 
thus been sought for almost a century. The value and usefulness 
of such measurements has been pointed out from time to time 
(e.g., von Hippel 1996) but conclusive or statistically significant 
observations are still lacking, which was a motivation for the 
present project. 

2. Radial velocities in open clusters 

The aim is to observe apparent radial velocities of main- 
sequence and giant stars in some stellar cluster(s) with sufficient 
precision to identify expected differences in gravitational red- 
shift between these groups. The ideal cluster should obviously 
be populous and have stars covering a large range of M/R ratios. 

Not all clusters are suitable. To assure good measuring preci- 
sion for spectroscopic radial velocity, the stars must have spec- 
tral lines that are both narrow and numerous, which excludes 
young clusters where main-sequence stars still are rapidly rotat- 
ing and produce broadened lines. Also, higher-temperature gi- 
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ants are unfit due to their few spectral lines and often complex 
atmospheric dynamics. Thus, the cluster must be sufficiently old, 
with main-sequence rotation slower than some 5-10 km s _1 , and 
must have a well developed red-giant branch with many cool gi- 
ants that are good targets for radial-velocity measurements (e.g., 
Setiawan et al. 2004). 

The relatively old cluster M67 (NGC 2682) at -850 pc dis- 
tance was chosen as a nearly ideal object. Besides having a well- 
developed red-giant branch, advantages include that its compo- 
sition is quite similar to solar (Tautvaisiene et al. 2000; Randich 
et al. 2006; Pace et al. 2008) and the age somewhat compara- 
ble: 2.6 Gy (WEBDA 2009). It has been well studied in various 
aspects (Pasquini et al. 2008; Sandquist 2004), including a cen- 
sus of binaries (Latham et al. 1992). Finally, we had already ob- 
served a homogeneous group of stars in M67 (Melo et al. 2001; 
hereafter Paper I), verifying the modest rotational velocities, and 
could use these data as a starting point. 

From previous radial-velocity data (e.g., Mathieu et al. 
1986), the internal velocity dispersion of the space motions 
within M67 was found to be m 0.5-0.8 km s _1 , similar to that 
deduced from proper motions (Girard et al. 1989). Other stud- 
ies suggest the expansion and rotation of the cluster core to be 
below 1 km s" 1 (Zhao et al. 1996; Loktin 2005). 

Since these numbers are comparable to the expected differ- 
ence in gravitational redshifts between dwarfs and giants, aver- 
ages over many stars will be required, and attention must be paid 
to calibrating radial-velocity measures between stars of differ- 
ent spectral types, identifying stellar multiplicity, possibly dif- 
ferent velocities for stars mass-segregated into different parts of 
the cluster, and the physics of line formation in dynamic stellar 
atmospheres. In any case, since the internal measuring precision 
of modern radial-velocity instruments is very much better than 
this dispersion, significant progress should be possible. 

3. Observations and data reduction 

An observing program was carried out with FEROS (Fiberfed 
Extended Range Optical Spectrograph; Kaufer et al. 1999) in 
three different runs. The first took place during FEROS com- 
missioning in October and November 1998, observing 28 main- 
sequence, turnoff, and giant stars. The radial velocities and 
Vsini were published in Paper I and are listed in the online 
Appendix for sake of completeness. Additionally 68, mostly 
main-sequence and turnoff stars, were observed over four nights 
during the first week of February 2002 (Run II) in parallel with 
other projects during Brazilian time. Finally in Run III, 66 addi- 
tional M67 members were observed in late February 2008. For 
the first two runs, the spectrograph was attached to the ESO 1 .5- 
m telescope, and during the third to the MPG/ESO 2.2-m tele- 
scope. In total 162 observations for 144 different stars were ob- 
tained. 

FEROS is a bench-mounted, thermally controlled, prism- 
cross-dispersed echelle spectrograph. It offers rather high res- 
olution (R ^ 48,000), and in a single recording with 39 echelle 
orders provides almost complete spectral coverage from ss 350- 
920 nm. On bright stars, FEROS has been demonstrated to be 
precise to better than 20 m s _1 (Setiawan et al. 2004). Our tar- 
gets are substantially fainter, but an exceedingly high precision 
is not required. The M67 velocity dispersion is on the order of 
1 km s~ 1 , so a precision of a few hundred meters in a single expo- 
sure is sufficient. Exposure times were set to have a S/N around 
30 or better. With this photon noise, the radial velocity (RV) pre- 
cision limit is of the order 20 m s _1 . The total error budget is 
discussed in Section 6.1. 



We aimed to observe many stars with different M/R ratios, 
to be able to precisely determine mean velocities for different 
groups of stars. Known binaries were excluded (M67 is rich in 
those), although some single-lined spectroscopic binaries likely 
are still present. We also observed some stars with deviant ve- 
locities which appear to be non-members. The list for the whole 
sample, including data from Paper I, is given in the Appendix, 
available online. 

4. Radial velocities and membership 

Radial velocities were computed as the Gaussian fit to the cross- 
correlation function (CCF) between the observed stellar spec- 
trum and a spectral template mask. This binary mask was based 
on a stellar spectrum, following the recipe in Baranne et al. 
(1979). As in Paper I, a K0 III synthetic mask (Baranne et al. 
1996) was used for all stars, assuring the same reference for all 
spectra. 

The mask is composed mostly of metallic lines, whose num- 
ber is sufficiently high (about 2000 lines in the 500-700 nm 
wavelength interval used) that we expect that line blends, which 
vary according to the stellar effective temperature and gravity, 
are substantially symmetric and do not change the average shift 
of the spectrum. The cross-correlation mask is built by select- 
ing the signal from spectral lines below a given residual flux. 
This cut is made in order to compromise between the number 
of points producing the cross-correlation signal and the width 
of the lines: too high a threshold (including points close to the 
continuum) might produce too wide CCFs, decreasing the mea- 
suring precision (see Baranne et al. 1996). For our mask a value 
of 0.7 was selected, i.e., only wavelength points with a residual 
intensity of less than 0.7 in the K0 III reference spectrum con- 
tribute to the digital mask. This implies that the selected lines 
have relatively large equivalent widths in the M67 giants. 

Only a few stars have repeated observations, and the com- 
parison between multiple measurements is shown in Figure [T] 
The average radial velocity of these stars in the first observing 
runs, made at the ESO 1.5-m, is 33.93 km s -1 , which compares 
very well with 33.96 km s _1 for the same objects observed at 
MPG/ESO 2.2-m, although the dispersion is quite large, around 
0.45 km s _1 . This number is, however, an upper limit to the 
measurement error, because among the 14 stars observed at the 
two telescopes, some probably are binaries. Three stars (S 1016, 
S1075, S1305) differ by more than 0.8 km s _1 , and once these 
are removed from the sample, the dispersion drops to 0.3 km s _1 . 
We consider this as a good estimate of the real uncertainty of a 
single measurement, though it might still be an upper limit. This 
value is similar to the dispersion in measurements of our RV 
standard star, HD 10700 (see section 5.1 and 6.1.2). 

Figure [2] shows the radial velocity distribution of all the ob- 
served stars. The M67 members are clearly visible, and their RV 
distribution can be well represented by a Gaussian distribution 
centered on RV=33.73, with cr=0.83 km s . We made a 2.5 
cr clipping and considered only the stars left as bona fide M67 
single members. Using this selection, a total of 33 stars were 
discarded. We note that S101 1, 1064, 1221, 1250, 1314, all clas- 
sified as binaries in Paper I, appear in the list of binary stars. Out 
of the Paper I and Latham et al. (1992) lists, only S986 was not 
recovered as a binary. The final sample consists of 1 10 stars and 
their color-magnitude diagram is in Figure [3] bona fide single 
members are marked as filled squares, stars discarded as empty 
squares. The targets are well distributed in the color-magnitude 
diagram, and numeorus evolved stars were observed to ensure 
good statistics also in that less populated group. 
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Fig. 1. Radial velocities for the 14 stars common to the differ- 
ent observing runs. Observations at the ESO 1.5-m telescope are 
black squares; those at the MPG/ESO 2.2-m are [red] triangles. 
The difference between the averages from the two telescopes is 
less than 30 m s , and the spread is 0.3 km s _1 after 3 suspected 
binaries are discarded. 



M67 is a well studied cluster, and among previous radial ve- 
locity studies, the most appropriate comparison is probably with 
the pioneering work by Mathieu et al. (1986), which is similar in 
size to ours (104 likely single members), with comparable (but 
slightly worse) RV precision. A difference is that their sample 
was limited to brighter magnitudes, since it included all proper 
motion candidates brighter than V=12.8, which corresponds to 
the region around the turnoff. Those authors used (as we do) a 
cross correlation technique, although their spectra were limited 
to a smaller spectral range, centered on 520 nm. We analyzed the 
Mathieu et al. (1986) radial velocities exactly in the same way as 
ours, finding a mean velocity 33.57 m s~'(cr= 0.78 km s~') and, 
anticipating the results of the next section, also their data do not 
show evidence of gravitational redshift. 



5. Analysis 

After cleaning the sample, we performed a gross split between 
'main sequence' and 'giants', defining as giants all stars along 
the subgiant and red-giant branches redder than B - V=0.7, and 
all the other stars as belonging to the main sequence. The gi- 
ants' radii are so large that gravitational redshift is expected to 
be negligible in these, and therefore their average RV should be 
bluer than the main sequence stars by several hundred m s -1 . 



Fig. 2. Radial velocity (RV) distribution of the observed stars 
in M67 (continuous line). A Gaussian distribution centered at 
Vr=33.73 and with cr=0.83 km s is superposed 



This simple split of the sample however does not produce the 
expected result, as evident from Figure |4j 

1. Contrary to what is expected from gravitational redshift pre- 
dictions, the two distributions are centered on almost the 
same velocity (33.75 km s -1 ) 

2. Main sequence stars and giants have substantially different 
velocity dispersions: ~ 0.90 and 0.68 km s _1 respectively. 

In proceeding with a more complete analysis, masses and 
radii were computed for all stars. While the mass range is small, 
the radii vary considerably. These were computed using stel- 
lar bolometric magnitude and effective temperature, calculated 
from the B - V vs. T e ff relation and the bolometric magnitudes of 
Alonso et al. (1996, 1999). A reddening of 0.041 (An et al. 2007; 
Taylor 2007) and a distance modulus of 9.61 (Pasquini et al. 
2008) for the cluster were assumed. For the masses, a fit between 
mass and luminosity was derived from evolutionary isochrones 
(Girardi et al. 2000) assuming solar metallicity and solar age. 
For the evolved stars, a flat initial mass of 1.32 M was adopted 
(however, small imprecisions in masses or radii do not influence 
our conclusions). 

The result is in Figure |5j RV as a function of the mass-to- 
radius {M/R) ratio. If gravitational redshift were the only mech- 
anism acting, one would expect a linear relationship, which is 
not present: the overall slope of the observed data can actually 
be fitted with a coefficient consistent with zero. 

A closer analysis suggests that a trend among the dwarfs 
(e.g., M/R >0.5) might be present, fitted by a RV = 0.99(M/R) 
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Fig. 3. M67 color-magnitude diagram with the observed targets. 
Empty boxes are stars discarded on the basis of their radial ve- 
locities differing more than 2.5 <x from the average, or known 
to be binaries. Filled squares denote stars retained as bona fide 
single members. 

law. However, given the small range in M/R (a factor 2 only) 
covered by the main-sequence stars, the uncertainty in this slope 
is large, and the slope is not statistically significant. Giants, if 
analyzed alone, show instead a negative slope of RV vs. M/R, 
also with a large error. 

5.1. Maximum likelihood estimation of the kinematic cluster 
parameters 

Following the approach laid out by Pryor & Meylan (1993) we 
performed a maximum likelihood estimation (MLE) of the kine- 
matic cluster parameters assuming a simple model including the 
effect of the gravitational redshift, and separating the two stellar 
groups of giants and dwarfs. We are interested in the cluster sys- 
temic velocity v c , the velocity dispersion of giants and dwarfs, 
cr g and cr d , respectively, as well as the dependence of the spec- 
troscopically measured velocity on (M/R) - presumably due to 
the gravitational redshift. We assumed that each individual mea- 
surement i was drawn from a probability density 



li(v c ,a,crj = 



1 



7 9 

oi + at 



exp - 



( v,- - v c - aPjY 
2(o- x 2 + a- 2 ) 



(1) 



Vj is the spectroscopically measured radial velocity, cr x the ve- 
locity dispersion of the group of stars the object i belongs to, cr m 
the measurement error on a single star for which we assumed it 
is the same for all stars and amounts to 0.3 km s -1 . For brevity, 



Fig. 4. Distribution of the measured radial velocities (RV) for the 
sample stars, sub-divided into 'giant' and 'dwarf sub-samples. 
Giants include all evolved stars with B-V>0.7. No difference is 
present in the central value of the two distributions. The dwarf 
distribution (dashed black line), shows a larger spread than the 
giants (904 vs. 680 ms" 1 ) 



we introduced P, = (M/R)j. a is a constant of proportionality 
which should come out as 0.636 km s -1 per unit M/R (in solar 
units) if the expected gravitational redshift was presenjM We as- 
sumed that there is no measurement uncertainty on the P,. This 
approximation is justified since the uncertainties are dominated 
by the spectroscopic ones which correspond to an uncertainty 
in Pi of ~ 0.5, much worse than the precision expected for Pj- 
values from isochrone fitting. 

We maximised the overall likelihood 

L(v c , a, cr g 2 , oj) = n fe{giants) Z;(v c , a, of) n je{dwarfs} / ; (v c , a, of), (2) 

i.e., we assumed that giants and dwarfs share the same systemic 
velocity, and exhibit an intrinsic line shift proportional to their 
mass to radius ratio. Giants and dwarfs can potentially have dif- 
ferent velocity dispersions - as is expected from the previous 
analysis. The ansatz for the likelihood above emphasises the role 
of intrinsic line shifts in the determination of the kinematic pa- 
rameters of the cluster, and the need for correcting for intrinsic 
shifts if one aims at a precise determination of the kinematic 
properties of cluster members. 

We followed the standard procedures for determining the 
most likely parameters and their covariance matrix, except that 
we used the variances instead of standard deviations as actual 



1 Here, we neglected the 3ms 1 redshift in the gravitational potential 
of the Sun at Earth orbit. 
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9 = (v c , a, eft , of). We obtained the best fitting parameters and 
their 1 o uncertainties as 



v c = 33.824 ±0.121 kms" 1 

a = -0.134 ±0.230 kms" 1 

cr g = 0.566 ±0.068 kms" 1 

o° d = 0.849 ±0.078 kms" 1 . 



(6) 



(M/R)/(M/R) 



The errors of the standard deviations were calculated from 
the errors on the variances applying the standard error propa- 
gation formula. From the covariance matrix we calculated the 
correlation coefficients among parameter pairs: -0.840 (v c , a); 
-0.067 (v c , of); 0.059 (v c , of); 0.093 (a, of); -0.081 (a, of); 
and -0.008 (of of). The only significant and actually strong 
correlation among the parameters is the one between v c and a. 
This is understandable since changes in the systemic velocity 
can to some extent be traded against changes of the intrinsic line 
shifts. The velocity dispersions in each group of stars are largely 
insensitive to the systemic velocity and the intrinsic line shifts. 
All in all, the MLE analysis confirms our findings that there is 
a clear difference (3.6 o) between the velocity dispersions of gi- 
ants and dwarfs, and a dependence of the intrinsic line shifts 
on (M/R) which deviates significantly (3.3 o) from the expected 
behaviour for a signal solely due to a gravitational redshift. 



Fig. 5. Examining the gravitational redshift predictions: Radial 
velocity vs. the mass/radius ratio (M/R). The black line is the fit 
to the data, while the dashed [red] line shows the linear depen- 
dence upon M/R expected from the gravitational redshift for- 
mula. Such a dependence is not seen, at least when all the stars 
are considered. Limiting the analysis to the 'main sequence' ob- 
jects (M/R > 0.5), a trend (though not statistically significant) 
might be present. 



6. Discussion 

Our data show two main results: that no obvious signature of 
gravitational redshift is seen in the data, and that the disper- 
sion of the radial velocities is higher in un-evolved stars than 
in evolved ones. We discuss these two results separately: 

6. 1 . Missing gravitational redshift 

As discussed above, a radial velocity shift of up to 600 m s is 
expected between main sequence stars and giants. Such a differ- 
ence is not present in our data as illustrated in Figures [4] and [5] 
Could we have missed the gravitational redshift because of some 
systematic effect? 



parameters in the minimisation process. This had mostly practi- 
cal reasons since some algebraic expressions became much more 
compact this way. For a single group of stars with velocity dis- 
persion cr x the best fitting parameters could be found analytically 

as 



a = 



(viPd-ivdiPi) 



{p'j) - (Pi) 2 

v c = (vd - a (Pi), 



and 



((Vi - v c - aPif) ■ 



(3) 



(4) 



(5) 



(.) denotes the (ensemble) average. We did not find an analytical 
solution for the case of two different velocity dispersions. For 
this case we solved the minimisation problem numerically by 
direct iteration on of and o^ as suggested by Pryor & Meylan 
(1993). With the best fitting parameters we evaluated the covari- 
ance matrix (-d 2 L/d8jd8jy l where 6 is the vector of parameters 



6.1.1. Template effects 

Given that the mask used is of spectral type K0 III, a better match 
between template and target star is expected for the sub-giants 
and giants than for main sequence stars. One clear effect of using 
the K0 III mask with G-type turn-off stars is that the resulting 
CCF is shallower than for the giants and sub-giants (Figure 3 
of Paper I). This is both because of the mismatch between star 
and template, and because of the generally weaker lines in dwarf 
stars; the depth of the CCF is thus related to the luminosity. 

Could the mismatch between luminosity class and mask have 
caused a velocity shift in the turnoff stars? To investigate this, 
we computed new cross-correlation functions to all stars using 
a G2 V mask instead. The mean RV difference was found to be 
a well defined value (30 m s -1 ) due to template zero-point mis- 
match. The RV difference distribution has an RMS of 120 m s _1 . 
A detailed inspection of spectra showing differences larger than 
100 m s _1 reveals that low S/N along with a not optimum cosmic 
ray cleaning are responsible. In addition, most of the discrepant 
spectra were recorded during the last season, when a simultane- 
ous Th-Ar lamp was used, contaminating faint spectra through 
its emission lines. If only data from the first two seasons are con- 
sidered, the RMS drops to 90 m s _1 , as in Fig. [6] 
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Fig. 6. Histogram of the differences of radial velocities obtained 
by using two different cross-correlation masks (K0 III giant and 
G2 V dwarf): RVk - RVc- The Gaussian fit gives a mean of 
-30 m s" 1 and an RMS of 90 m s -1 . 



While these tests show some limits in the precision, the use 
of a different mask also shows that there is no offset introduced 
by the choice of the mask differentially affecting various groups 
of stars. However, an extra error of about 100 m s _1 due to poor 
cosmic -ray cleaning and ThAr contamination has to be added to 
the formal error of the cross-correlation fit. 

We also tested different mask characteristics, such as depen- 
dence of the measured RV upon the wavelength range, finding 
that values retrieved from the 500-700 nm mask are within a 
few tens of m s _1 from those found for the 'full' 400-1000 nm 
one, for different types of stars. We therefore conclude that the 
mask and wavelength range used do not induce systematic ef- 
fects which could cancel the gravitational redshift signal. 

6.1.2. Spectrograph drifts 

Variations in temperature and atmospheric pressure during the 
observing night can cause spectral drifts of a few hundred m s _I 
(e.g., Platais et al. 2007). This effect is seen in our data collected 
in the third season, where the simultaneous calibration fiber was 
used. Unfortunately, the first two datasets did not make use of the 
simultaneous calibration fiber and the spectrograph drifts could 
not be removed. 

In order to estimate the amplitude of the drifts caused by 
changes in atmospheric conditions we used 52 observations of 
r Cet (HD 10700) taken from the SACY survey (Torres et 
al. 2006) and from our own observations. These spectra were 



H 10 13 14 16 

V Magnitude 



Fig. 7. Residuals of measured radial velocity VR as function of 
stellar apparent magnitude. The residuals show a typical increase 
with magnitude as expected from increasing noise for fainter 
stars. The effect might also be due to a dynamically hotter pop- 
ulation of dwarfs, or to some observational bias. 

reduced and analyzed in a similar way to the present M67 
stars. They cover a time span of 3000 days from Oct. 1999 to 
Feb. 2008. The mean RV is -16.37 km s 1 with an RMS of 
0.3 km s _1 . The best fit, excluding 6 extreme values, gives a 
similar mean, and a sigma of 0.18 km s _1 . Thus we believe that 
the drift error is about 200 m s . 

Simulations of the effect of photon noise in the RV computa- 
tion show that already with a S/N of 10-30, radial velocity errors 
are as low as a few tens of m s _1 (Paper I; Bouchy et al. 2001). 
However those simulations only consider white (photon) noise, 
and (as discussed above) there are additional sources of errors 
associated with low S/N. In Figure|7]we plot the scatter diagram 
as a function of stellar magnitude and indeed the behavior is that 
expected, if the noise was not negligible. 

We investigate if one can reproduce the missing gravitational 
trend and the scatter seen in Figure[7]accounting for the different 
sources of spread in the VR. First, the intrisic velocity dispersion 
of the cluster of 0.8 km s _1 was considered. A Gaussian drift er- 
ror drawn from a distribution with an RMS of 200 m s and 
a photon noise of 150 m s for a V=14 star were added. This 
noise was scaled as the square root of the flux for a given mag- 
nitude (i.e., l/SN, assuming that photon noise is the limiting 
factor for faint targets). 

We then computed the difference of the center of the RV dis- 
tribution of main sequence and giant stars: this is always below 
100 m s _1 . Moreover, the RMS of both distributions is clearly 
dominated by the cluster velocity dispersion. An additional grav- 
itational shift of 600 m s" 1 added to the distribution of sub-giant 
and giant stars can be retrieved in all cases in spite of their lim- 
ited number. 

We cannot so far identify any instrumental or observational 
effect which could cancel out the expected gravitational redshift 
signal, and we conclude that our observations show the absence 
of evidence for a gravitational redshift trend. 

Nordstrom et al. (1997) briefly discuss a similar case, be- 
cause they found a difference of 0.4 km s _1 (dwarfs redder than 
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Fig. 8. Synthetic reproduction of the observed residual RV vs. 
magnitude diagram of Figure [7] As input for the points, our var- 
ious sources of noise were used, together with the velocity dis- 
persion inside the cluster. 



giants) between the average RV of giants and main sequence 
stars in the open cluster NGC 3680, almost exactly what is pre- 
dicted from gravitational redshift. NGC 3680 is an almost com- 
pletely evaporated cluster, with only a few bona-fide members 
left. The paucity of the sample precludes certain conclusions, but 
could suggest that our results about M67 may not be generalized 
to other clusters. On the other hand, the fact that no gravitational 
redshift is found by analyzing the published data by Mathieu et 
al. (1986) shows that this result is not limited to our observations 
only, but it is a characteristic of this cluster. 



6.1.3. Cluster dynamics 

Even if open clusters are not known to rotate or to have com- 
plex dynamical topologies, we shall analyze whether the missing 
gravitational signature could be due to a combination of complex 
cluster dynamics and bias selection, namely that the observed gi- 
ants and dwarfs are situated in different spatial regions of the 
cluster, perhaps affected by different global velocities. In this 
case the stars would not share the same systemic velocity and 
this bias could conspire to cancel the gravitational redshift sig- 
nature. For M67, there is no report of complex dynamics, with 
the exception of a claim by Loktin et al. (2005) of a possible 
expansion of the cluster core. No evidence of rotation is found 
in the analysis of the Mathieu et al. (1986) data by the same au- 
thors. In order to investigate effects related to cluster dynamics 
we have analyzed the velocity distribution of the observed stars 
as a function of R (distance from the cluster center) and 6 (an- 
gular position) for giants and dwarfs, finding no evidence of RV 
variations with these parameters (Figure EJ. We note in passing 
that we do not find evidence for a larger dispersion at large radii. 
We therefore conclude that complex cluster dynamics coupled to 
a selection bias cannot explain the non-detection of gravitational 
redshift. 



Fig. 9. Angular and radial distribution of measured radial veloc- 
ities (black squares denote giants; [blue] dots are main sequence 
stars). No evidence of complex dynamics or dynamical heating 
at larger distances is observed. 



6.1.4. Convective shifts 

Dynamic effects in stellar atmospheres may produce net shifts in 
photospheric absorption lines, and we cannot exclude that lines 
in main-sequence stars could be more blueshifted than in giants, 
thus (partially) canceling the expected gravitational redshift sig- 
nature. 

Such a differential shift need not be present in all spectral 
lines, but only (in a statistical sense) in the deeper ones. In fact, 
for a given metallicity, typical lines (metallic neutral ions) are 
generally deeper in cool giants than in dwarfs. By using the same 
mask for all stars, we use lines with different equivalent widths in 
the two classes of stars, and relatively deeper in the giant spectra. 

Since single spectra of each star have a rather low S/N ratio, 
we have constructed two 'average' low-noise spectra: one from 
64 selected dwarfs (called 'dwarf') and one for 20 selected gi- 
ants ('giant'). The spectra were first normalized in flux at 600 
nm, and then averaged by calculating the medians. Given the 
low velocity dispersion of the cluster stars, these co-added spec- 
tra are not appreciably broadened. The S/N ratio of these spectra 
is about 200 in the region of interest for the 'dwarf spectrum 
and about 140 for the 'giant' one. The co-added spectra still 
show - within the uncertainties - the same radial velocity, for 
'giant' and 'dwarf alike. We selected a number of lines used 
for chemical analysis of solar metallicity stars (Pasquini et al. 
2004), inspecting them visually in the co-added spectra. Out of 
the original list, 52 lines passed this quality check, and were fit- 
ted with a Gaussian function, using the ARES package (Sousa 
et al. 2008). We first tested with the IRAF task Splot the con- 
sistency between the automatic program and the manual fitting. 
The central wavelengths in the 'dwarf and 'giant' spectra were 
recorded, and the DV = V g i ant - Vd wa rf is plotted in Figure [10] 
as function of line equivalent widths. For 'giant' lines weaker 
than about 4 pm (40 mA), the measures were not very reliable, 
and such weak lines are not used. On the other hand, lines with 
equivalent widths above 15 pm (150 mA) were discarded to min- 
imize possible contamination by blends. In both plots of Figure 
10 a trend is present: the stronger the line, the redder the shift 
in the giant spectrum. The observed relative shifts are up to 
500 m s _1 , comparable to the expected gravitational redshift. 
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Fig. 10. Difference between measured radial velocities of se- 
lected lines in the median 'giant' (black squares) and 'dwarf 
spectra ([red] triangles). A trend of ARV vs. equivalent width 
[mA] is present: stronger lines show smaller (negative) differ- 
ences between dwarfs and giants. For strong lines, the expected 
gravitational redshift probably is canceled by atmospheric dy- 
namical effects. 



We looked for, but did not find any dependence with other quan- 
tities (e.g., wavelength). 

Formal fits to the trends give: DV=2.1*EWd„arfs-350 ms" 1 
and DV=5.0*EW g i ants -418 ms -1 . The trends are only marginally 
significant (correlation coefficients 0.25 and 0.35). The formal 
error in measurements of single lines is about 50 ms~', and we 
believe that the observed spread is real. The lines used are mostly 
from neutral elements (dominated by Fe) with different charac- 
teristics, but also a few Fe II. With higher S/N spectra for single 
stars, it would be possible to further investigate trends between 
different chemical elements and line parameters, but we believe 
that for our current data this would be an over-interpretation. 

Wavelength shifts induced by convective motions have in 
particular been studied for lines in the Sun and in Procyon, 
where it has been possible both to measure and to model intrin- 
sic lineshapes under very high spectral resolution (Dravins 1982; 
Allende Prieto et al. 2002a); a successively increased blueshift 
for weaker lines is a signature of the three-dimensional photo- 
spheric hydrodynamics (e.g., Asplund et al. 2000). Also, in other 
cool stars, Allende Prieto et al. (2002b) and Ramirez et al. (2008) 
observed how lines from neutral species get more blueshifted as 
they become fainter. 

Since we are able to compare dwarfs and giants of the same 
cluster, some information can be added: the trend seems to differ 



between main sequence and evolved stars. In the M67 evolved 
stars, the blueshift is less than in turnoff ones, and the expected 
gravitational redshift is seen only for a few of the faintest lines. 
The blueshift variation with equivalent width is larger in giants 
than in dwarfs, as suggested by the trend in Figure 10 A caveat, 



however, is that further parameters may enter, in particular line 
broadening by stellar rotation which differs among dwarfs and 
giants: rotational broadening changes line asymmetries and may 
affect the velocity signal measured by cross correlation (Dravins 
& Nordlund 1990). 

This discussion suggests a plausible explanation for our non- 
detection: gravitational redshift is largely canceled because our 
correlation mask is dominated by rather deep lines, which are 
more blueshifted in main sequence stars than in giants because 
of photospheric convective motions. 

6.2. Convective line shifts in 3D model atmospheres 

In order to investigate the apparent absence of gravitational red- 
shift, synthetic spectral lines were computed from 3-dimensional 
and time-dependent model atmospheres. We synthesized 15 ar- 
tificial Fe I lines, spanning the range of excitation potential and 
equivalent width of the observed ones. All lines were placed at a 
wavelength of 620 nm, which is the average of the observed Fe I 
lines in the range 586-672.7 nm. The 3D model atmospheres 
providing the thermal and kinematic information for the spec- 
tral synthesis were taken from the CFIST 3D model atmosphere 
grid (Ludwig et al. 2009). For representing the 'dwarf and 'gi- 
ant', we chose for each a single 3D model atmosphere. Table [2] 
summarizes their basic properties. 

The wavelength shifts of the synthesized lines were mea- 
sured applying the same procedure as for the observed ones, 
i.e., by fitting a Gaussian profile. In contrast to the observations, 
the artificial lines were available at very high spectral resolution 
(s; 2 x 10 6 ) and signal-to-noise ratio so that their shifts with re- 
spect to the laboratory wavelength could be calculated with high 
fidelity. For the 32 observed Fe I lines we calculated convective 
line shifts by interpolating in the grid of synthetic lines accord- 
ing to the observed equivalent width and excitation potential of 
the line. The typical statistical uncertainty to which the wave- 
length position of an individual line is predicted is (from pre- 
vious experience) 30ms~'. This means that for the difference 
between wavelength shift of the same line between two models 
we expect a typical precision of 42 ms -1 . The statistical uncer- 
tainty is related to issues like how well the acoustic oscillations 
present in the computational box are averaged out. The example 
illustrates that we further expect that the line shifts among lines 
of the same model are correlated, hence we do not expect that the 
uncertainty of the predicted shifts among giant and dwarf drops 
as the square-root of the number of lines used to characterise it. 
Thus, we expect a modeling uncertainty of about ±40 ms -1 for 
the line shift difference between the giant and dwarf - not a large 
value but also not completely insignificant. 



Figure 11 shows the predicted convective shifts in dwarfs 



and giants. Before determining the line shift, the synthetic spec- 
tral line profiles were broadened with a Gaussian profile of 
<x 3.9kms _1 which corresponds to the spectrograph resolu- 
tion, combined with the macroturbulence-like velocity disper- 
sion stemming from the motions of individual cluster stars. In 
addition, it was assumed that the dwarf rotates with a projected 
rotational velocity of 5kms _1 while the giant exihibits no ro- 
tation. We find for the difference of the shifts of lines between 
the giant and dwarf a minimum of 62 a maximum of 211 and 
an average of 146 ms -1 . We note that the correlation between 
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Table 2. 3-dimensional radiation-hydrodynamics model atmospheres. 



Model 


^eff 


logg 


[Fe/H] 
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Vrms 
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[K] 


[cm/s 2 ] 




[Mm] 
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M 


% 
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dwarf 


5 782 
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2.27 
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12 
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giant 
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2.50 
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160 
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5 
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18.9 


3.04 
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Notes. "Model" is the model's name used in this paper, r er y the effective temperature, logg the gravitational acceleration, [Fe/H] the metallicity, 
l x = l y the linear horizontal size of the square-shaped computational box, Z z its vertical extent, N x = N y the number of grid points in the horizontal 
directions, N z the number in vertical direction, T the duration of the simulated time series, N„b ra the number of equivalent frequency points 
considered in the solution of the radiative transfer equation, SI^/I the relative spatial white light intensity contrast at stellar disk center, v rms the 
temporal and spatial average of the vertical velocity at Rosseland optical depth unity, "Modelcode" is an internal identifier of the model sequence. 
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Fig. 11. Calculated convective wavelength shifts for the 32 stud- 
ied Fe I lines for the dwarf (crosses) and giant (squares) models. 
Negative velocities correspond to absolute blueshifts (gravita- 
tional redshifts are left out). 



the strengths of lines in the giant and dwarf is not very tight. 
However, as seen in Fig. [TT[ we find a trend of the convective 
shift with equivalent width which is continuous between dwarfs 
and giants: on average the stronger Fe I lines are less blueshifted 
in the giant. 

Figure [12] depicts the calculated differences of the shift of 
the same line between dwarf and giant, comparing them to ob- 
servations. For the synthetic lines, assumptions about spectro- 
graph resolution and stellar rotation were the same as before. 
An average gravitational redshift was now included, so that the 
giant synthetic points were blueshifted by 382 ms -1 (applying 
an average M/R of 0.73 for dwarfs and 0.13 for giants, as for 



the stars used to compute the average spectra). The agreement 
with the observations is now quite good, with a difference of 
only ~ 140 ms -1 between the synthetic and the observed aver- 
age shifts. 

However, the observed spectra show effects not predicted by 
the modeling: dependence of the shift on the equivalent width, 
and a considerable spread in the measurements. Nevertheless, 
we conclude that convective line shifts are the most likely ex- 
planation as to why the expected gravitational redshift signature 
was not observed. 

Among possible causes of remaining discrepancies between 
models and observations, we note that the chosen 3D models 
may not be the optimal representatives of the 'dwarf and 'giant' 
spectra. Their choice was dictated by availability. For example, 
the iron abundance determined from the dwarf model is approx- 
imately solar but for the giant it turns out to be 0.5 dex greater. 
This may indicate that the giant model is hotter than the typical 
observed giant. For the dwarfs, solar temperature is appropriate 
for our fainter stars, while higher temperatures are expected for 
the stars closer to the turnoff. Indeed a similar Fe abundance is 
obtained by adopting 6000 and 4600 K for dwarfs and giants 
respectively (using ID models). However, we do not expect the 
line shifts to sensitively depend on stellar parameters (as long 
as the line strength is matched) so we argue that our approach 
should give a reasonably accurate estimate of the expected line 
shift differences. 

Another possibility is the presence of noise and biases in our 
combined spectra; in order to test that we repeated the analy- 
sis using the Kitt Peak solar flux spectrum (Kurucz et al. 1984) 
and the Arcturus spectrum (Bagnulo et al. 2003). Although here 
we can only test the slope and the spread around the relation- 
ship, this gives qualitatively similar results to M67, although the 
slope is shallower (~ 2 m s~7 m A instead of 5) and the scatter is 
smaller (100 m s instead of 240 m s ). 

Being the first attempt of such a test, we believe it is en- 
couraging. In presence of higher S/N data it will be possible 
to test other effects, such as the dependence of the line shifts 
on other parameters (element, excitation potential, ionization, 
wavelength). One could also apply the same technique to clus- 
ters of different metallicity and chemical composition, and to 
work out comparison grids for theoretical models. 

6.3. Dynamically hot dwarfs 

The observation that dwarfs are dynamically hotter than giants 
might seem counter-intuitive, because giants are known to be 
intrinsic RV variables, with variations of up to a few hundreds 
m s _] (e.g., Setiawan et al. 2004), and their intrinsic RV variabil- 
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Fig. 12. Comparison of observations to synthetic wavelength 
shifts of Fel lines, computed from 3-D hydrodynamic mod- 
els for dwarfs and giants. The computed differences between 
the lineshifts in giant and dwarf models is shown as crosses, 
and the observations as squares. The calculated shifts include 
the average expected contribution of the gravitational redshift of 
382 ms- 1 



ity will add to their RV dispersion. In addition, main sequence 
stars have lower masses, and in case of mass segregation the gi- 
ants should be (as they are in M67) mostly concentrated in the 
middle of the cluster, near the center of the potential well. On 
the other hand, in an equipartition regime, giants should be dy- 
namically cooler, because they are more massive. 

The proper-motion dispersion of the bright M67 stars as 
measured by Girard et al. (1989) is almost 0.2 mas/yr. At a dis- 
tance of 836 pc (Pasquini et al. 2008) it corresponds to a radial 
velocity dispersion of 0.8 km s _1 , comparing very well with the 
dispersion of RV=0.83 km s measured from all stars and also 
with the 0.68 km s _I dispersion of the giants only. The compar- 
ison of the radial velocity distributions for both populations in 
Figure|4]suggests that the dwarfs are dynamically hotter than the 
giants in M67, i.e., they have a larger velocity dispersion. 

In order to check whether this might originate by some 
observational bias, we performed the following test: we again 
simulated the radial velocity distribution for the dwarfs using 
the same error parameters and a population of 100 stars as in 
Figure [8] The RMS of this artificial distribution was then com- 
puted, and the procedure repeated 5000 times. The result is a 
distribution of RMSes, centered at the velocity dispersion of the 
cluster (800 m s -1 ) with values between 0.6 and 1.0 km s _1 . If 
in the simulations we now turn off the spectrometer drift error 



Fig. 13. Distribution of velocity dispersions for 5000 realizations 
of simulated populations of 1000 (dark) and 100 stars (light), 
using actual parameters for observational errors. 



(Sect. 6.1.2), while maintaining all other parameters, we end up 
with similar figures. This already indicates that the observed ve- 
locity dispersion does not depend on measurement errors, but 
rather on the number of objects. 

As a further test, we set the drift error to its initial value and 
generated the RMS distributions for two populations: one with 
N=100 members (corresponding to our real case) and another 
with 1000 members. The result is in Figure 13 With 1000 mem- 
bers, the velocity dispersion is well defined with a center value 
of 0.8 km s _1 and 0.73 to 0.88 km s peak-to-peak. Therefore, 
the observation that the M67 dwarfs have a larger velocity dis- 
persion than the giants seems not caused by our measurement 
errors but possibly is an effect of low number statistics. 

We note that Nordstrom et al. (1997) found a similar re- 
sult for NGC 3680, with the dwarfs showing a dispersion of 
0.67 km s compared to 0.36 of the giants. Also, analysis of 
the Mathieu et al. (1986) data (which, however, contain a limited 
number of main sequence stars) is in perfect agreement with our 



findings, giving cr(RV)- 0.89 and 0.66 km s for dwarfs and 
giants respectively. 

Since we cannot exclude that the effect is real, it is worth 
mentioning that a similar result was found by Zhao et al. (1996) 
and Loktin et al. (2005) from analyzing proper motion data. 
Zhao et al. (1996), analyzing a sample of about 200 stars, found 
that the dispersion of proper motions of M67 bright stars (V< 
13.5 ) is significantly smaller than that of fainter objects (13.5 
< V < 14.5), by a factor ~1.3. They also found that the brighter 
(more massive ) stars are more concentrated towards the cen- 
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ter of the cluster. Also Loktin et al. (2005) found a similar re- 
sult, with a much larger sample. The larger sample allows them 
to subdivide the color-magnitude diagram into some sub-sample 
with giants, and another with turnoff and upper main sequence 
stars. These two subsamples closely correspond to our division 
between giants and dwarfs. For dwarfs they find dispersion ve- 
locities on average almost 1 .3 times higher than for giants (their 
Table 3). 

We also note that the tangential and radial velocity disper- 
sions derived by Loktin et al. (2005) are substantially larger than 
our radial- velocity dispersions: 0.84 and 0.77 km s _1 for the 
giants; 1.01 and 1.04 km s _1 for dwarfs ('radial' in their work 
refers to polar coordinates in the sky plane). Loktin et al. (2005) 
adopt a distance to the cluster of 904 parsecs. Using our pre- 
ferred distance of 831 pc, the radial and tangential dispersions 
become (0.77, 0.71) and (0.93, 0.96) km s" 1 , respectively for 
giants and dwarfs, fully compatible with the observed radial ve- 
locity dispersions of 0.68 (giants) and 0.9 km s -1 (dwarfs). This 
result indirectly confirms the shorter distance scale for this clus- 
ter that emerges from recent works (cf. Pasquini et al. 2008). 

It is remarkable that in both our data and those by Mathieu 
et al. (1986), the ratio between the radial velocity dispersions of 
main sequence and evolved stars is 1.3, the same as found by 
Zhao et al. (1996) and Lotkin et al. (2005) for the dispersion be- 
tween fainter and brighter stars. Zhao et al. (1996) interpret these 
results claiming that the dynamical energy of the more massive 
stars has been passed in the time evolution of the cluster to the 
less massive objects, which are now escaping the cluster. 

While a full dynamical analysis of the cluster is beyond the 
scope of this work, the possibility of joining precise radial ve- 
locity and proper motion data seems attractive for detailed future 
investigations. 

M67 is by now the best studied open cluster when combin- 
ing both, the number of stars observed and the precision of the 
measurements. Substantially larger samples of radial velocity 
measurements have been gathered in globular clusters (see e.g. 
Sommariva et al. 2009), but they are dynamically much more 
complex objects and with a larger velocity dispersion. 

7. Conclusions 

A search was carried out to find the expected differences in grav- 
itational redshift between the spectra of dwarf and giant stars in 
the M67 open cluster, whose stars should share the same sys- 
temic velocity, and where the average radial motion for differ- 
ent types of stars thus should produce the same Doppler shift. 
However, despite adequate measuring precision, the expected 
signature was not seen. Investigations of the reasons for ab- 
sence of signal - in particular following high-resolution spectral 
synthesis from 3-dimensional hydrodynamic model atmospheres 
- indicate that likely reasons are different amounts of convec- 
tive blushifts in different types of stars. An enhanced convective 
blueshift in, particularly, those stronger spectral lines in dwarf 
stars that carry a larger weight in the cross-correlation techniques 
used for radial-velocity measurement, appear adequate to cancel 
the expected gravitational redshift. 

This study thus also illustrates the limitations of cross- 
correlation techniques for high-accuracy lineshift measure- 
ments. To understand the origin of wavelength shifts that 
are only a small fraction of any spectral-line width, a cross- 
correlation measure over an extended spectral range may not 
be adequate, but rather one needs to examine (groups of) indi- 
vidually selected line profiles measured with sufficient spectral 
fidelity to reveal also their intrinsic asymmetries and amount of 



rotational broadening. These can then be compared to synthetic 
spectral lines from 3-D atmospheres complete with convective 
blueshifts and gravitational redshifts. Not many current spec- 
trometers have adequate performance for such tasks, but work 
is in progress towards realizing such high-fidelity spectroscopy 
with an accurate wavelength calibration at the largest telescopes 
(Dravins 2010; Pepe et al. 2010; Pasquini et al. 2010). With these 
instruments it will be possible to extend the present study to less 
explored stars, for instance in the metal-poor and in the metal- 
rich regimes. 
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